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Abstract 

Using 58 million J /ip and 14 million 'ip{2S) events collected by the BESll 
detector at the BEPC, branching fractions or upper limits for the decays 
J/ijj and ip{2S) AAvr^ and AAr] are measured, and the decays of J/iJj and 
ipi^S) to nKgA + c.c. are observed and measured for the first time. Finally, R 
measurement data taken with the BESII detector at center-of-mass energies 
between 3.7 and 5.0 GeV are fitted to determine resonance parameters of the 
high mass charmonium states, ip{3770), '?/'(4040), ip{4160), and '?/'(4415). 

The Beijing Electron Collider is being upgraded to a two-ring collider 
(BEPCII) with a design luminosity of 1 x lO^^cm"^ s~^ at 3.89 GeV and will 
operate between 2 and 4.2 GeV in the center of mass. With this luminosity, 
the new BESIII detector will be able to collect, for example, 10 billion J/ip 
events in one year of running. BEPCII and BESIII are currently nearing 
completion, and commissioning of both is expected to begin in mid-2008. 

1 Introduction 

In this paper, some recent BESII results are reported based on 58 million J/ ip 
and 14 million ip{2S) events collected by the BESII detector at the BEPC, 
and the status of BESIII/BEPCII is summarized. For much more detail, see 
the references. 



2 Recent Results 

2.1 J/V^ and V^(25) decays into AAtt^ and AA?7 

The isospin violating decay J/ip — >■ AAtt^ was studied by DM2 [1] and 
BESI [2], and its average branching fraction is B{J/ip AAvr^) = (2.2 ± 
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Table 1: Measured branching fractions or upper limits at the 90% confidence 
level (C.L.). 



Channels 


Number of events 


MC efficiency (%) 


Branching fraction (x 10 ^) 


J/V' ^ AAttO 


< 11.2 


0.75 


< 0.64 


J/'f/i — > AA?7 


44± 10 


1.8 


2.62 ±0.60 ±0.44 


ip{2S) AAir" 


< 7.0 


2.5 


< 0.49 


^{28) ^ AAt] 


< 7.6 


2.9 


< 1.2 


J/%1> -> E+TT-A 


335 ± 22 


2.3 


7.70 ±0.51 ±0.83 




254 ± 19 


1.8 


7.47 ±0.56 ±0.76 



0.6) X 10 ^ [3]. However, the isospin conserving process J/ip ^ AAt] has not 
been reported, and there are no measurements for AAvr^ and AA?7 decays of 

Here J/ip AAvr^, J/ip AAr], ^(2S) ^ AAtt^, and ^^(2^) ^ AArj, 
where A decays to vr^p and tt" and r] to 77, are studied. Candidate events 
must have four good charged tracks and at least two photons, two protons 
identified using particle identification, a satisfactory four constraint kine- 
matic fit, and a np mass consistent with the A mass. For J/i(j —>■ AAvr^, the 
decay lengths of A and A in the x — y plane must be larger than 0.05 m. 

Possible backgrounds are studied using MC simulation. The most serious 
one is from J/iJj —>■ S°7r°A + c.c, which contains AAvr^ with an additional 
photon in the final state. The branching fraction for this decay has not 
been previously measured. Since direct measurement of J/ip ^ S°7r°A + c.c. 
is difficult, we measure the branching fractions of its isospin partner and 
estimate the branching fraction assuming isospin symmetry. 

The histogram in Fig. [1] shows normalized backgrounds from all back- 
ground channels, and the dashed line in the figure shows the 7r° signal from 
MC simulated J/ijj AAvr^. The dat a in Fig. [T] are consistent with back- 
ground, and the upper limit on the number of 7r° events from J/ip ^ AAvr^ 
is determined. Figure [2] shows the the 77 invariant mass distribution for 
J/ip ^ AAr] candidates, where a clear rj signal is observed. Figure [3] shows 
the 77 invariant mass distribution for 'ip{2S) AAtt" and AAr/ using similar 
selection criteria, and no significant tc^ or rj signals are seen. 

Table [T] lists the results for J/tp and ip{2S) decay into AAvr^ and AArj, 
as well as J/tp S+vr^A + c.c. Except for J/ip AAvr'' and J/ip ^ 
E+7r~A + c.c, the results are first measurements. Interestingly, the result 
of J/ip AAvr" presented here is much smaller than those of DM2 and 
BEST [1,2]. Previously, the large contaminations from J/ip TPti^A + c.c. 
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Figure 1: Invariant mass distribution 
of M(77) for J/%1) AA7r° AA77 
candidates (dots with error bars) and 
normahzed backgrounds (sohd his- 
togram). The dashed curves shows 
the 7r° signal from MC simulated 
J/i/j AAttO. 
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Figure 2: Fit to the 77 invariant 
mass distribution of J/ip ^ AAt] 
AA77 candidate events. Dots with 
error bars are data, the hatched 
histogram is the normalized back- 
ground, and the solid histogram is 
the fit to data using a histogram of 
the signal shape from MC simulation 
plus a second order polynomial for 
background. 



and J/tp S+TT^A + c.c. were not considered, resulting in a large J/ip ^ 
AAtt^ branching fraction. The small branching fraction of J/ip ^ AAvr*^ and 
relatively large branching fraction of J/ip AArj measured here indicate 
that the isospin violating decay in J/ decays is suppressed while the isospin 
conserving decays is favored, which is consistent with expectation. For more 
detail, see Ref. [4]. 

2.2 J/i/^ and '0(25) decaying to nKgA + c.c. 

In 2004, BESII reported the observation of an enhancement X(2075) near 
the threshold of the invariant mass spectrum of pA in J/ip pK^A decays. 
The mass, width, and product branching fraction of this enhancement are 
M = 2075 ± 12 (stat.) ± 5 (syst.) MeV/c^, F = 90 ± 35 (stat.) ± 9 (syst.) 
MeV/c2, and B{J/^ K~X)B{X pA+c.c.) = (5.9±1.4±2.0) x lO'^Jb], 
respectively. The study of the isospin conjugate channel J/ip —>■ nK^A is 
therefore important not only in exploring new decay modes of J/ip but also 
in understanding the X(2075). 

J/ip and ip{2S) — > nK^A with Kg tt+tt" and A pir^ (and c.c.) final 
states contain four charged tracks and an undetected neutron or anti-neutron. 
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Figure 3: The 77 invariant mass 
distribution for candidate ip{2S) — > 
77AA events. Dots with error bars 
are data, and the histograms are MC 
simulated signal events. The arrows 
indicate the vr^ and 77 signal regions. 



Figure 4: The invariant mass spec- 
tra of (a) Air^, (b) nK'l, and (c) 
A?2(Arr), as well as (d) the Dalitz 
plot for candidate J/ip ^ nK^K + 
c.c. events after all selection criteria. 
The crosses show the sideband back- 
grounds. 



We require the candidate events to have four charged tracks with total charge 
zero. Secondary vertex fitting is used to identify the tt+tt" from the and 
the pTT from the A, and their masses are required to be consistent with those 
of the parent particles. To reject backgrounds from channels without a K'g 
or A, we require L^.y(A), the distance from the reconstructed A vertex to the 
event origin, to be larger than 5 mm and Lxy{K^) > 5 mm. To suppress 
background and improve the resolution, a one constraint (IC) kinematic fit 
with a missing neutron is applied under the J/ip ^ pnn^n^n^ hypothesis, 
and Xic < 5 is required. 

The invariant mass spectra of Ai^"^, nKg, and An(An), as well as the 
Dalitz plot for all selection requirements are shown in Fig. |H In the AKg 
invariant mass spectrum, an enhancement near AKg threshold is evident, as 
is found in the AK mass spectrum in J/ip ^ pK~A [6]. The X(2075) signal 
which was seen in the pA invariant mass spectrum in J/ip ^ pK~~ A is not 
significant here. Taking into account the systematic error, the upper limit of 
the near-threshold enhancement X(2075) in the nA threshold is B{J/iIj 
irgX(2075)) ■ 5(X(2075) ^ nA) < 4.9 x IQ-^ (90% C.L.). Considering 
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the isospin factor, the branching fraction upper hmit for B{J/%1) K^X) ■ 
B{X nA + c.c.) is not inconsistent with that for B{J/ip KX) ■ B{X — > 
pA + c.c.) [5]. 

An N* state at around 1.9 GcV/c^ in the Ai^'° invariant mass spectrum 
and A* states at around 1.5 and 1.7 GeV/c^ in the nK^ invariant mass 
spectrum are present. A larger data sample and a partial wave analysis are 
needed to analyze these states. 

We use the same criteria to select ■0(25') — > nK^h. + c.c. events from the 
BESII sample of 14M ip{2S) events. The branching ratios obtained are: 

Br{J/ij nK%K + c.c) = (6.46 ± 0.20 ± 1.07) x 10"^ 
5r( J/V' ^ nK%K) = (3.09 ± 0.14 ± 0.58) x 10"^ 
Br{J/i) nK%k) = (3.39 ± 0.15 ± 0.48) x 10"^ 
Br{'il){2S) nK%K + c.c) = (0.81 ± 0.11 ± 0.14) x 10"^ 

The ratio of the branching ratios of ip{2S) and J/^ decaying to nKgh. + c.c, 
gft=(12.6±3.5)%, obeys the "12%" rule [7]. For more detail, see Ref. [8]. 



2.3 V?(3770), i/?(4040), i/?(4160) and V?(4415) resonance 
parameters 

The total cross section for hadron production in e+e" annihilation is usually 
parameterized in terms of the ratio i?, which is defined as i? = (7(e+e~ — > 
hadrons)/(7(e"'"e~ — fi'^fi~), where the denominator is the lowest-order QED 
cross section, a{e^e~ A*"*"/^") = 0"°^ = 47r«^/3s. At the open flavor 
thresholds where resonance structures show up, R measurements are used to 
determine resonance parameters. For the high mass charmonium resonances, 
the V'(3770) was measured by MARK-I [9], DELCO [10], MARK-II [11] and 
BES [12] [13]; the ^'(4040) and V(4160) were measured by DASP [14]; and 
the V(4415) was measured by DASP [14] and MARK-1 [15]. 

The most recent and precise R measurements between 2-5 GcV were 
made by BESll [16]. Experimentally, R for both the continuum and the 
wide resonance region is given by 

R = ^fegrf ~ ^bg 

<^'^t^Letrgehad{i + Sobs) ' 

where N^^d number of observed hadronic events, N^g is the number 

of the residual background events, L is the integrated luminosity, (1 -|- Sobs) 
is the effective correction factor of the initial state radiation (ISR) [17] [18], 
^had is the detection efficiency for hadronic events determined by Monte Carlo 
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Table 2: The resonance parameters of the high mass charmonia in this work 
together with the values in PDG2004 [19], PDG2006 [3] and K. Seth's eval- 
uations [20] based on Crystal Ball and BES data. 







i/;(3770) 


^-(4040) 


7/;(4160) 




M 

(MeV/c2) 


PDG2004 
PDG2006 
Ur> (betnj 
BES (Seth) 
BES (this work) 


3769.9±2.5 
3771.1±2.4 

3771.4±1.8 


4040±10 
4039±1.0 
4Uo i ±Z 
4040±1 
4038.5±4.6 


4159±20 
4153±3 

4i0i±4 

4155±5 
4191.6±6.0 


4415±6 
4421±4 

44Z0±D 

4455±6 
4415.2±7.5 


Ttot 
(MeV) 


PDG2004 
PDG2006 
CB (Seth) 
BES (Seth) 
BES (this work) 


23.6±2.7 
23.0±2.7 

25.4±6.5 


52±10 
SOiblO 
85±10 
89±6 
81.2±14.4 


78±20 
103±8 
107±10 
107±16 
72.7±15.1 


43±15 
62±20 
119±16 
118±35 
73.3±21.2 


r 

^ ee 

(keV) 


PDG2004 
PDG2006 
CB (Seth) 
BES (Seth) 
BES (this work) 


0.26±0.04 
0.24±0.03 

0.18±0.04 


0.75±0.15 
0.86±0.08 
0.88±0.11 
0.91±0.13 
0.81±0.20 


0.77±0.23 
0.83±0.07 
0.83±0.08 
0.84±0.13 
0.50±0.27 


0.47±0.10 
0.58±0.07 
0.72±0.11 
0.64±0.23 
0.37±0.14 


(5 (degree) 


BES (this work) 





133±68 


301±61 


246±86 
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Figure 5: The fit to the R values in the high mass charmonium region. The 
dots with error bars are the updated R values. The solid curve shows the 
best fit, and the other curves show the contributions from each resonance 
Rbw, the interference Rint, the summation of the four resonances Rres = 
Rbw + Rint, and the continuum background Rcon respectively. 
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simulation without bremsstrahlung being simulated, and etrg is the trigger 
efficiency. 

In the previous analysis [16], the determination of R values was done us- 
ing PDG04 [19] resonance parameters for the high mass resonances. Here, 
the analysis uses the data to determine the resonance parameters. The de- 
termination of R values and resonance parameters are intertwined; the factor 
(1 + Sobs) in Eq- dl]) contains contributions from the resonances and depends 
on the resonance parameters. Therefore, the procedure to calculate {l + 5obs) 
requires a number of iterations before stable results can be obtained. We 
perform a global fit over the entire center-of-mass energy region from 3.7 
to 5.0 GeV covering the four resonances, '?/'(3770), ■?/'(4040), ipi^AlQQ) and 
'?/'(4415), and include interference effects among the resonances. We also 
adopt energy-dependent full widths, and introduce relative phases between 
the resonances. 

The resonant parameters of the high mass charmonia determined in this 
work, together with those in PDG2004, PDG2006 and the results given in 
Ref. [20] are listed in Table El The updated R values between 3.7 and 5.0 
GeV and the fit curves are illustrated in Fig. [5l 

It is worth noting that the change of the resonance parameters affects 
the effective initial state radiative correction factors, and thus affects the R 
values. In general the relative difference is within 3%, and for a few energy 
points the maximum difference is about 6%. Our resonance parameter results 
are in agreement with the previous experiments in most cases, but large 
differences are observed in some of the parameters, such as the mass of the 
'i/'(4160). This is mainly due to the reconsideration of the radiative correction 
factors, and the inclusion of interferences between the resonances. This work 
is preliminary; for more detail, see Ref. [21]. 

3 BEPCII and BESIII 

In 2003, the Chinese Government approved the upgrade of the BEPC to 
a two-ring collider (BEPCII) with a design luminosity approximately 100 
times higher than that of the BEPC. This will allow unprecedented physics 
opportunities in this energy region and contribute to precision flavor physics. 

3.1 BEPCII 

BEPCII is a two-ring e^e~ collider that will run in the tau-charm energy 
region {Ecm = 2.0 — 4.2 GeV) with a design luminosity of 1 x 10^^ cm~^s~^ at 
a beam energy of 1.89 GeV, an improvement of a factor of 100 with respect 
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Figure 6: Schematic view of the BESIII detector. 



to the BEPC. This is accomplished by using multi-bunches and micro-beta. 
The upgrade uses the existing tunnel. 

The 2024 meter long linac has been upgraded with new klystrons, a new 
electron gun, and a new positron source to increase its energy and beam 
current; it can accelerate electrons and positrons up to 1.89 GeV with an 
positron injection rate of 50 mA/min. Its installation was completed in the 
summer of 2005. 

There are two storage rings with lengths of 237.5 meters. The collider 
has new super-conducting RF cavities, power supplies, and control; super- 
conducting quadrupole magnets; beam pipes; magnets and power supplies; 
kickers; beam instrumentation; vacuum systems; and control. The old dipoles 
are modified and used in the outer ring. Electrons and positrons will collide 
at the interaction point with a horizontal crossing angle of 11 mrad and bunch 
spacing of 8 ns. Each ring has 93 bunches with a beam current of 9170 mA. 
The machine is already providing a high flux of synchrotron radiation at a 
beam energy of 2.5 GeV. 



The BESIII detector consists of a berylium beam pipe, a helium-based small- 
celled drift chamber, Time-Of-Flight counters for particle identification, a 
CsI(Tl) crystal calorimeter, a super-conducting solenoidal magnet with a 
field of 1 Tesla, and a muon identifier using the magnet yoke interleaved 
with Resistive Plate Counters (RPC). Fig. [H] shows the schematic view of the 
BESIII detector, including both the barrel and endcap portions. 



3.2 BESIII 
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3.3 Physics in the tau-charm energy region 

The tau-charm energy region makes available a wide variety of interesting 
physics. Data can be taken at the J/ip, ip{2S), and ip{3770), at r threshold, 
and at an energy to allow production of Dg pairs, as well as for an R-scan. 



Table 3: Number of events expected for one year of running. 



Physics 


Center-of-mass 


Peak 


Physics 


Number of 


channel 


Energy 


Luminosity 


cross 


Events per 




(GeV) 


(10^3 cm-2 s"^) 


section (nb) 


year 




3.097 


0.6 


~ 3400 


10 X 10^ 


T 


3.67 


1.0 


~ 2.4 


12 X lO'' 




3.686 


1.0 


~ 640 


3.0 X 10^ 


D 


3.770 


1.0 


~ 5 


25 X 10*^ 


Ds 


4.030 


0.6 


- 0.32 


1.0 X 10*^ 


Ds 


4.140 


0.6 


~ 0.67 


2.0 X 10*^ 



BEPCII and BESIII are in the final stage of assembly, and commissioning 
will begin in summer 2008. The design luminosity of BESIII is 1 x 10^3 cm~^. 
Clearly BESIII with higher luminosity will contribute greatly to precision 
flavor physics: Vcd and Vcs will be measured with a statistical accuracy of 
1.6%. D^D^ mixing will be studied and CP violation will be searched for. 
Table [3] gives the numbers of events expected during one year of running at 
various energies. Huge J/ip and ipi^S) samples will be obtained. The rjc, 
Xcjy and he can be studied with high statistics, and the pvr puzzle will be 
studied with better accuracy. For more detail, see Refs [22,23]. 
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